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Abstract 

We study the stabihty of some strongly frustrated antiferromagnetic spin lattices in high 
magnetic fields against lattice distortions. In particular, we consider a spin-s anisotropic 
Heisenberg antiferromagnet on the square- kagome and kagome lattices. The independent 
localized magnons embedded in a ferromagnetic environment, which are the ground state 
at the saturation field, imply lattice instabilities for appropriate lattice distortions fitting to 
the structure of the localized magnons. We discuss in detail the scenario of this spin-Peierls 
instability in high magnetic fields which essentially depends on the values of the exchange 
interaction anisotropy A and spin s. 
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1 Introduction. Localized magnons 

The square-kagome and kagome Heisenberg antiferromagnets have attracted much interest in 
recent times because of interplay between quantum fluctuations and frustrated lattice structure 
^miSlSIS]- Recently, for a wide class of frustrated spin lattices, which includes both the square- 
kagome and kagome lattices, exact eigenstates consisting of independent localized magnons in a 
ferromagnetic environment have been found |nil3IH|- They become ground states if the saturating 
magnetic field is applied and lead to a macroscopic jump in the zero-temperature magnetization 



Figure 1: Square-kagome lattice with one distorted square (top) and kagome lattice with one 
distorted hexagon (bottom) which can host localized magnons. The parts of the lattices before 
distortions are shown by dashed lines. All bonds in the lattice before distortions have the same 
length. 

curve just below saturation. Moreover, very recently we have examined a field-tuned instability 
of the square-kagome and kagome spin lattices with respect to lattice distortions through a 
magnetoelastic mechanism [S] reporting rigorous analytical results completed by large-scale exact 
diagonalization data for lattices up to = 54 sites TD . That study was restricted mainly to 
the spin-^ isotropic Heisenberg antiferromagnet. Since often one meets kagome materials with 
s > ^ (for instance, spin-| kagome-like compound Ba2Sn2ZnCr7pGaio-7p022 with a comparably 
small exchange constant of about 37 ... 40 K ) it is useful to go beyond the case s = ^ and 
A = 1. In this paper we extend our previous results for the spin quantum numbers s > ^ and 
exchange interaction anisotropy A ^ 1 demonstrating that a scenario of the discussed earlier 
magnetic-field induced lattice instability may essentially depend on the values of s and A. 
To be specific, we consider the square-kagome lattice (Fig. ^ top) and the kagome lattice 



(Fig. ^ bottom). The Hamiltonian of quantum spins reads 

H=Y,Jnm {sts- + S-S+) + A<S^ J - hS\ (1) 
(nm) ^ ^ 

Here the sum runs over the bonds (edges) which connect the sites (vertices) occupied by spins for 
the lattice under consideration, Jnm > are the antiferromagnetic exchange constants between 
the sites n and m, A > is the anisotropy parameter, h is the external magnetic field, and 
= Y^j^ is the z-component of the total spin. We assume that all bonds in the lattice 
without distortion have the same length and hence all exchange constants have the same value 
J. 

From Refs. [HI El El E] we know that independent localized one-magnon states embedded in 
a ferromagnetic background are exact ground states of the Hamiltonian at saturation field 
for the considered models. 

More specifically, by direct computation one can check that 

|1) = ^7 d'^ ~ 1,S,'S,S) - |S,S - 1,S,S) 

+ \s, s,s - l,s) - \s,s,s,s - l)) \ . . . s . . .) (2) 

and 

|1) = (|s — 1, s, s, s, s, s) — \s,s — 1, s, s, s, s) + \s, s,s — 1, s, s, s) 
— |s, s,s,s — 1, s, s) + \s, s, s, s, s — 1, s) — \s, s, s, s, s, s — 1)) \ . . . s . . .) (3) 

are one-magnon eigenstates of the Hamiltonian on the square-kagome and kagome lattices. 
Here s or s — 1 denote the value of s^, the omitted site indices in the first multiplier in the r.h.s. 
of Eq. ((21) (Eq. Q) run along a square (hexagon) trapping cell, and the second multiplier in 
the r.h.s. of Eqs. and ((SJ, | . . . s . . .), stands for the embedding fully polarized ferromagnetic 
environment. The corresponding energies {h = 0) of the one-magnon states ((2) and ((SJ are 

El = -2s J + 2s(2s - 1)A J + 2s(4s - 1)A J + {2N - 12)s^A J (4) 

and 

El = -2s J + 2s(3s - 1)AJ + 2s(6s - 1)AJ + {2N - 18)s'^AJ. (5) 

We separate in Eqs. @, © the contributions to the energy from those bonds which form 
a magnon trapping cell (first and second terms), from the bonds connecting this cell with 
the environment (third terms) and from the ferromagnetic environment (fourth terms). The 
considered lattices may contain n = 1, . . . ,nmax localized iiV or iiV 



for the square-kagome and kagome lattices, respectively. Each magnon decreases the total 
by one and a localized magnon state with n independent magnons has 5^ = sN — n. In the 
presence of an external field h ^ the energy E{S^, h) can be obtained from the energy without 
field E{S^) through the relation E{S^, h) = E{S^) - hS\ 

Under quite general assumptions it was proved [01 El that the localized magnon states have 
lowest energies in the corresponding sectors of total S^. As a result, these states become ground 
states at the saturation field. More specifically, the ground-state energy in the presence of a 
field is given by Eo{h) = E^i^iiS^) — hS^ and the ground-state magnetization 5^ is determined 
from the equation h = E^i^{S^) — £'min('S'^ — !)• Since for = sN, . . . , sN — nmax the localized 
magnon states are the lowest states, one has 

E„,in{S') = 2Ns'^AJ - nei = 2Ns'^AJ - Nsei + eiS' (6) 

where 

ei = 2s(l - (2s - 1)A) J - 2s(4s - 1)AJ + 12s^AJ = 2s(l + 2A) J (7) 

and 

ei = 2s(l - (3s - 1)A) J - 2s(6s - 1)AJ + ISs^AJ = 2s(l + 2A) J (8) 

for the square-kagome and kagome lattices, respectively. Due to the linear relation between 
i^min and © one has a complete degeneracy of all localized magnon states at the saturation 
field h = hi = ei, i.e. the energy is 2Ns'^AJ — Nsei at h = hi for all sN — "/imax ^ < sN. 
Consequently, the zero-temperature magnetization jumps between the saturation value sN 
and the value sN — {sN — ^N) for the square-kagome (kagome) lattice. 

Finally, we mention that the localization of magnon in a finite area has some relation to the 
flat one-magnon dispersion. Some rigorous results for the flat-band electronic systems have been 
reported in Refs. [HI [HI HSj . 

2 Spin-Peierls instability 

To check the lattice stability of the considered systems with respect to a spin-Peierls mechanism 
we assume a small lattice deformation which preserves the symmetry of the cell which hosts 
the localized magnon (in this case the independent localized magnon states remain the exact 
eigenstates) and analyze the change in the total energy. The corresponding deformations are 
shown in Fig. ^ For the square-kagome lattice (Fig. ^ top) the deformations lead to the 
following changes in the exchange interactions: J — > (1 + \/^6)J (along the edges of the square) 



and J —^ {1 — i-v/2(\/3 — l)S)J (along the two edges of the triangles attached to the square), 
where the quantity 6 is proportional to the displacement of the atoms and the change in the 
exchange integrals due lattice distortions is taken into account in first order in 5. For the 
kagome lattice (Fig. ^ bottom) one has J — > (1 + 5) J (along the edges of the hexagon) and 
J ^ (1 — ^S)J (along the two edges of the triangles attached to the hexagon). Note, that the 
lattice distortions shown in Fig. ^ correspond to 5 > 0. The magnetic energies (01) and Q are 
changed by distortions by (^-2V2s + 2^/2s{2s - 1)A - ^^s(4s - 1)A^ 6 J and -s(2 + A)6J, 
respectively. Note that the linear with respect to 5 change of the magnetic ground-state energies 
has a simple reason: the magnetic ground-state energy for sN — n-max ^ < sN contains a 
sum of contributions coming from local lattice regions each of which varies linearly with 6. The 
elastic energy increases in harmonic approximation by 2 (6 — \/3) (square-kagome) and by 
9a6^ (kagome). The parameter a is proportional to the elastic constant of the lattice. The 
changes of total energies due to distortions read 

AEi = ^-2^25 + 2V2s{2s - 1)A - ^^ ^(4^ - 1)A^ 5J + 2(6- Vs^ a5^ (9) 

and 

AEi = -s(2 + A)5J + 9a<5^ (10) 

where for n independent localized magnons trapped by distorted cells these results have to be 
multiplied by n. Here comes the first important conclusion. For the kagome lattice AEi UlUj) 
decreases at small > for any A > and any s, thus implying a spin-Peierls instability with 
the distorted trapping cells like in Fig. ^ bottom. In contrast, for the square-kagome lattice with 
s > ^ the sign of 5 which provides a decrease of AEi for the small lattice distortions depends 
on the anisotropy parameter A: 5 > if A < A*(s) (with A*(l) 2.21748, A*(|) « 0.92171 
etc; A*(s) ~ ^ as s ^ 00) but 5 < if A > A''(s). This can be understood result of 
interplay between the contributions of transverse and zz correlations along different bonds for 
the square-kagome lattice geometry. 

Let us discuss the scenario of spin-Peierls instability for the square-kagome lattice with 
A < A*(s) and for the kagome lattice. Consider a magnetic field above the saturation field hi. 
For the corresponding fully polarized ferromagnetic state a lattice distortion is not favorable. 
Decreasing h till hi the homogeneous ferromagnetic state transforms into the "distorted magnon 
crystal". Further, numerical analysis suggests that both spin-1 systems exhibit a magnetization 
plateau between hi and /i2 < hi at = sN — rimax- Calculating the plateau width A/i = /ii — /12 
for spin-1 isotropic (i.e. A = 1) finite systems of = 24, 30 (square-kagome) and N = 27, 36, 45 



(kagome) for the undistorted lattice (/i2 is obtained by h2 = -Emin('S'^ = sN — n^^y^ — Ey^\^{S^ = 
sN 

~ f^max ~ 1)) and using a jj finite-size extrapolation we find some indications for a finite 
A/i ~ 0.5J (A/i ~ O.IJ) for the square-kagome (kagome) lattice in the thermodynamic limit. 
We also mention here that on the basis of general arguments jl61 I17j one may expect the 
magnon crystal state to have gapped excitations that is related to the magnetization plateau at 

S'-" = sN - n-max- 

Now the question arises whether the lattice distortion under consideration is stable below 
this plateau, i.e., for < sN — rimax- We discuss the question again for spin-1 finite systems 
of size N = 24 (square-kagome) and = 36 (kagome) with rimax distorted squares/hexagons. 
We calculate the magnetic energy for zero and small distortion parameter 5 for different values 
of . Adopting for the magnetic energy the ansatz 

i?min(5^ 6) = i?min(S^ 0) + Ad^ (11) 

and taking S of the order of 10~^ we estimate the exponent p from the numerical results. 
Evidently, the lattice may become unstable if the magnetic energy 1)11(1 decreases with the 
exponent p < 2 whereas p > 2 indicates lattice stability. Interestingly, in the sector of 
just below the magnon crystal, = sN — rimax — 1 (i-e. as h becomes smaller than /12), the 
spin-1 square-kagome and kagome lattices show different behavior. For the kagome lattices with 
= 36 we do not find a lattice instability. Note that this is in agreement with the situation 
for s = ^ ^Oj. Contrary to that, for the square-kagome lattice with = 24 we find a lattice 
instability for 6 > also in the sector = sN — nmax — 1 = 19 in case of small enough 
anisotropy A < A(l) « 1.6. Again if A > A(l) the favorable lattice distortion is characterized 
by (5 < 0. Since we know from Ref. [TOl that for the spin s = ^ square-kagome lattice the 
spin-Peierls instability can be observed even for lower < sN — n^s.x — 1; we check this also for 
the spin s = 1 square-kagome lattice with N = 24. As for = 19 we find a lattice instability 
for 6 > only for small enough anisotropy A < 1.0 for = 18 and A < 0.51 for = 17, 
i.e. anisotropy parameter A below which this instability occurs becomes smaller when (i.e. 
the magnetic field h) is diminished. Thus, we arrive at the second important conclusion. For 
the spin s = 1 kagome lattice with any A > the spin-Peierls instability is favorable only for 
sN 

"'^max ^ 5*^ < sN and the distortion disappears for /i < /12. In contrast, the distorted spin-1 
square-kagome lattice with the sufficiently small anisotropy parameter < A < A(l) remains 
stable for smaller h < h2- The spin-1 square-kagome lattice with A(l) < A < A*(l) exhibits 
more intricate behavior: the parameter 5 characterizing the lattice distortion at the saturation 
changes its sign for h < h2. 



As in the spin-^ case JOl the saturation field in the distorted lattice is shifted to higher 
values that provides a hysteresis phenomenon in the vicinity of saturation field. 

3 Concluding remarks 

There is an increasing number of synthesized quantum kagome magnets often with large values 
of spin (unfortunately, however, the available at present materials are not perfect kagome Heisen- 
berg antiferromagnets) and we may expect that further such materials will be synthesized. With 
our analysis we are pointing out that the efforts in this direction are worthwhile also because of 
a new effect: the spin-Peierls instability in a strong magnetic field. From the point of view of 
possible experiments, it is important to bear in mind the following remarks. First, large values 
of s leads to decrease of a relative plateau width For instance, we have ~ 0.11 for s = ^ 
[ini and ^ 0.08 for s = 1 (undistorted square-kagome, A = 1); ^ w 0.023 for s = ^ [101 
and ^ 0.017 for s = 1 (undistorted ka gome, A = 1). Thus, we need materials with not 
too large s for which the plateau width is not extremely small. Second, we need materials with 
comparably small exchange constant J to have experimentally accessible saturation field hi. 
Third, the presented consideration refers to ideal kagome geometry. We have not discussed the 
effects of deviation from the perfect kagome geometry on the predicted spin-Peierls instability 
in high magnetic fields. Nevertheless we may expect that the effects of localized magnon states 
will survive (see, for example, Ref. ^Hl)- Fourth, the reported analysis is performed within the 
frames of the adiabatic treatment of the spin-Peierls instability. The problem becomes much 
more difficult when the nonadiabatic effects are taken into account. 

To summarize, we have examined a spin-Peierls instability in strong magnetic fields for 
two frustrated spin-s anisotropic Heisenberg antiferromagnets hosting independent localized 
magnons demonstrating cumulative effects of exchange interaction anisotropy, spin value and 
lattice geometry. 
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